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SUMMARY 

In this work, we investigate the communication reliability for diffusion-based molecular communication, 
using the indicator of bit error rate (BER). A molecular classified model is established to divide molecules 
into three parts, which are the signal, inter-symbol interference (ISI) and noise. We expand each part 
separately using molecular absorbing probability, and connect them by a traditional-like formula. Based on 
the classified model, we do a theoretical analysis to prove the feasibility of improving the BER performance. 
Accordingly, an adaptive threshold variation (ATV) algorithm is designed in demodulation to implement the 
goal, which makes the receiver adapt the channel condition properly through learning process. Moreover, the 
complexity of ATV is calculated and its performance in various noisy channel is discussed. An expression 
of Signal to Interference plus Noise Ratio (SINR) is defined to verify the system performance. We test some 
important parameters of the channel model, as well as the ATV algorithm in the simulation section. The 
results have shown the performance gain of the proposal. Copyright © 2010 John Wiley & Sons, Ltd. 
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1. INTRODUCTION 

Molecular communication (MC) is an attractive domain which uses the molecules as the carriers to 
permit the information exchange between nano-devices at nanoscale [1]. MC has the advantage of 
energy-saving and no-radiation, which is more suited in biological environment compared with the 
traditional communication. Rapid growth of nano-technology promotes the manufacture of bio¬ 
inspired nano-devices, which are able to perform basic tasks, including sensing, actuating and 
computing. Limited by size and power, multiple nano-devices are essential to work together as 
a network to execute a complex task, composing the internet of nano-things. One typical case 
is the body area network, in which the nano-devices enjoy a bright prospect on nano-healthcare 
[2], and will gradually take place of traditional applications, such as the wireless sensor network 
[3,4]. Some another potential application domains include environmental monitoring, industry and 
military equipment, respectively described in [5]-[7]. 

There are some significant differences between traditional communication and MC, such as the 
information carrier, transmission speed, range, noise source and so on. MC could also be divided 
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into various forms according to specific rules. Comparing with other forms in MC, diffusion-based 
MC is inspired by drift flows or relies on molecular thermodynamics movement, which derives 
from the natural cases and is supported by fluid and molecular mechanics. Diffusion-based MC is 
worthy to be investigated not only in communication, because it could be applied in many useful 
applications. One typical example is the drug delivery [8,9], in which particles are released from the 
sickness, diffused in the blood medium for accurate positioning. Hence, we adopt diffusion-based 
MC as the research scenario in our paper. 

Communication reliability is an important issue all the time. Considering the randomness and 
uncertainty of molecular movement, MC suffers a more serious problem compared with traditional 
communication. At present, bit error rate (BER) is chosen as the indicator of the reliability in 
MC generally. For example, [10] does a brief investigation of reliability using BER in diffusion- 
based MC, but lacks deep analysis as the specific channel characteristics. [11] proposes the forward 
error correction codes over MC and verifies by BER. [12] designs the receiver with different signal 
detection technologies, in which the performances of those proposals are shown using BER. 

In our work, on-off key (OOK) modulation is adopted as most works do. Assuming transceivers 
are strictly time synchronous, one time slot is utilized to transmit one bit. We first introduce the basic 
system model and communication process. Then, based on the molecular absorbing probability, we 
propose a classified model to divide the molecules into signal, inter-symbol interference (ISI) and 
noise branches. Inspired by traditional communications, we establish a brief formula to connect the 
three branches. We expand each branches in detail, and give the expression of BER. 

In our work, to improve the performance of BER, we propose an algorithm named adaptive 
threshold variation (ATV). A relevant work is introduced in [13], where an adaptive transmission 
rate method is designed by fixing the receiver threshold and altering the transmission rate of 
transmitter. Inversely, our ATV algorithm fixes the transmission rate and alters the receiver threshold 
adaptively. We first prove the feasibility of the ATV design through 3 Propositions in theoretical 
analysis. Then we design ATV in detail. In ATV, the threshold varies based on the knowledge of the 
previous bits, which is related on the channel condition. The threshold variation is similar with the 
learning process, that make the receiver adapt the channel for the demodulation. We calculate the 
complexity of the ATV and discuss its operation under various noisy conditions. In the simulation 
section, we deduce the expression of signal to interference plus noise ratio (SINR). The results 
show that the optimal threshold is not at the general midpoint of the molecule quantity per bit, 
which is consistent with the theoretical analysis. ATV algorithm is verified in the SINR-BER curve, 
indicating that there is a significant performance improvement. In addition, we show the threshold 
variation for various algorithm parameter and under different communication conditions, to give the 
guidance of algorithm setting. 

The remainder of this paper is organized as follows. In section 2, we list the related work in 
this field. In section 3, we introduce the basic system model of the diffusion-based MC. Then, 
the molecular classified model is proposed in section 4. In this section, signal, ISI and noise are 
respectively analyzed. In section 5, the ATV algorithm is analyzed, designed and discussed in detail. 
The numerical results are presented in section 6, where we evaluate the BER performance in the 
analytical model, as well as the algorithm performance. We conclude this paper in section 7. 


2. REEATEDWORK 

Molecular communication is promising to be utilized in body environment, implemented based 
on the biological channels. According to different mediums, MC proceeds with different forms. 
Some potential forms under study includes diffusion-based MC [14], neural communication [15], 
molecule motor MC [16], blood vessel communication [17], bacteria-based communication [18], 
etc. Among those forms, diffusion-based MC is a typical one, which is similar with the traditional 
wireless communication. Major efforts in the area of diffusion-based MC focus on the study of 
physical layer at present. 

Diffusion-based molecular communication is usually a composite physical process, driven by 
two different motivations. For the former, molecules in the environment perform the Brownian 
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Figure 1. The diffusion-based MC system 


movement [14]. For the latter, molecules are driven by active drift flow [19]. Most of the works 
adopt synchronous mechanism, i.e., time slots are used to transmit bits. Differently, [20] proposes 
an asynchronous method for MC. For channel estimation, some works focus on the molecular 
concentration in the environment, which follows the well-known Fick’s second law [1]. Another 
works investigate the random walk of solo molecule, and model the distribution by winner process 
[19]. To evaluate the communication performance, some important indicators of diffusion-based 
MC system, like channel capacity [21], delay [22], bit error rate [10], are exploited based on the 
knowledge of traditional communication and information theory. 

Diffusion-based MC suffers serious interference problem in time-slotted system because of the 
randomness of molecule moving, even the regularities of distribution could be tracked. Inter-symbol 
interference (ISl), one of the typical interference in MC, are caused if molecules of other bits impact 
the current bit [23]. Some methods to eliminate ISl are designed accordingly. For instance, [24] 
tries to use the enzyme to remove the redundant molecules, and alleviate the interference between 
different slots. A decision back decoding method is proposed in [25], which could eliminate ISl 
based on the calculation of previous bits. The noise of MC is determined by several factors. [26] 
does a deep analysis of noise for diffusion-based MC. The choice of carrier molecules is a major 
one, because there may exist various background molecules in the environment, such as the human 
body. A specific molecular type should be chosen as the information carrier to alleviate the effect 
of those environmental molecules. Another main factor is interference from another nano-devices, 
for a nanoscale device to device (D2D) communication, this could be classified into noise source in 
MC. 

Many modulation methods of MC are extended from the traditional ones, including AM [27], FM 
[28], PM modulation [20], etc. In MC, AM modulation is based on the concentration of molecules, 
which maps potential of traditional communication. ”01” bits are obtained via comparison between 
concentration and a decision threshold. Among AM modulations, OOK is the simplest way and 
adopted in most MC literatures. FM modulation utilizes different types of molecules, which map 
various frequency in traditional communication. The interference between types of molecules 
comes from the possible chemical reactions, so it’s important to choose the appropriate types. PM 
modulation of MC is based on the release time of molecules. In one time slot, the molecules could 
be released at different time, that lead the difference of the phase. 


3. SYSTEM MODEL 

In this paper, one transmitter (TX) and one receiver (RX) nano-devices are considered, with the 
fixed distance denoted by r, as Eigure 1 shown. The essential communication process contains 
modulation, emitting, propagation, receiving and demodulation. We assume that the transceiver is 
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Figure 2. Transceiver slots mapping diagram 


Strictly time synchronous. A major goal of MC is to realize the reliable transmission of binary 
strings, which is denoted by h{i) in our model. Each bit is emitted in a transmitting slot nt and is 
respected to be absorbed in a corresponding receiving slot rir. Those time slots have the fixed length, 
given by r. One example diagram of transeeiver slots mapping is shown in Figure 2. 


3.1. Modulation and emitting 

Modulation is the first proeess of the eommunication. On one hand, restrained by the tiny size and 
limited power of nano-devices, it’s not easy to achieve some complex eommunieation technologies. 
One the other hand, many natural cells usually release specifie molecules when communication is 
necessary, and not release molecules during the rest time. So the on-off key (OOK) modulation 
eould fit many natural eases. We express the OOK modulation as. 


Ntx{i) 


M b(i) = ”l” 
0 b(i) = ”0” 


( 1 ) 


Where Ntx{i) is the released molecule quantity at the beginning of the slot. M moleeules are 
emitted for bit ”1”, no molecules are emitted for bit ”0”. Note that molecule quantity could not be a 
negative value, so it is a unipolar modulation. 


3.2. Propagation 

The movement of the moleeules in the environment eontains two folds. One is driven by active flows, 
the other is the random brownian movement, of which the composite mobility eould be modeled 
by winner process. From the macroscopical view, molecular concentration is a useful indicator 
in MC, similar with the potential in traditional eommunieation. Molecular concentration depends 
on several faetors, including time, spatial position, transmission rate, channel condition, etc. In 
one-dimensional space, the molecular concentration distribution follows the Fick’s second law [1], 
expressed as. 


= £)Ac(a;,f)-f r(0,t) (2) 

Where c(x, t) is the molecule concentration at location x and time t. Ac(x, t) means the sum of 
the second derivatives of c{x,t). r{0,t) is the transmission rate of the molecules at transmitter, of 
which the discrete form is given by r(0, it) = Ntx{i)- D is the diffusion coefficient related to some 
environment factors, given by. 


D = 


ksT 

6'!rr]RH 


( 3 ) 
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Where ks is the Boltzman constant, T is the environment temperature, rj is the dynamic viscosity 
of the fluid, and Rh is the hydraulic radius of the molecules. 

For various r(0,f), different solutions are available for (2). The OOK modulation is adopted 
as above said, the solution of (2) could be the integral of molecular transmission rate and Green 
function g{x, t) [29], which is a function of time and position, i.e., 

c{x,t) = r{0,t) * g{x,t) (4) 

Note that g{x, t) could also be regarded as the position distribution of a single molecule. In one¬ 
dimensional space, it is given by. 


g{x,t) 


(47rDf)3/2 AD? 


(5) 


3.3. Receiving and demodulation 

Let Vr be the maximum absorbing space of the receiver. Apart from the noise, we calculate the 
molecules in this space during slot as, 

p{i+l)T 

Nrx(.i)= / c{r,t)*VRdt (6) 

J IT 

In our model, we consider a realistic species of biological receiver, named ligand-based receiver. 
The specific ligand receptors are distributed uniformly on the surface of the receiver. It works in 
case that the carrier molecules are matched with the receptors. In other word, the receptors and 
signal molecules shall be chosen particularly, to be matched biologically. There are many cases 
in nature. For example, in neuron network, the receiver neural cells use specific ligand-based 
chemoreceptors to receive neurotransmitter in gap junction [30]. The receptors bind and release 
molecules continually, and (6) is rewritten as. 


Nr 


S) = [ 

J in 


(z+1)t 


c(r, t)*VR*a*Q 


dt 


(7) 


Where a and b are the binding and releasing rate of receptors. Q is the density of receptors on the 
surface of the receiver. 

The receiver demodulates the information according to the absorbed molecule quantity Nrx{i) of 
(7). Like traditional AM demodulation, a threshold is needed, that is fixed in general. We denote it 
by Nr. Bit ”1” is obtained if Nrx{i) is greater than threshold, else bit ”0” is obtained. Hence, we 
get the error probability of bit for bit ”1” as. 


N'T — 1 

Te(*) = ^ p{Nrx{i) = k) ( 8 ) 

/t-0 

Similarly, the error probability of bit for bit ”0” is, 

oo 

Pe(*) = X] P(^rx{i) = k) (9) 

k=NT 

In (8) and (9), p],{i) andpg(j) respectively means the error probability transmitting of ”1” and ”0”. 
They correspond with the sum of probability which the estimation is ”0” and ”1”. The expression 
ofBER is. 


Pe = p{b{i) = 0)p°(i) +p{h{i) = l)pl{i) (10) 

Where p{h{i) = 0) and p{b(i) = 1) are the probability of transmitting ”0” and ”1”. In many cases, 
we assume that they are equal as 0.5. 
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4. MOLECULAR CLASSIFIED MODEL 

In this section, we first give the absorbing probability of single molecule. Then, a classified model of 
molecules for demodulation is proposed, to divide the received molecules into signal, inter-symbol 
interference and noise parts. We establish a traditional-like formula to connect them, and expand 
them in detail. 

4.1. Molecular absorbing probability 

Considering that the distance between the transceiver is short (usually several /rm), we assume that 
the molecules live long enough in the environment before absorbed by the receiver. The distribution 
of a single molecule is given by (5). We calculate its accumulative density function G{x, t) as. 


G{x, t) 


g{x, t)dt 



( 11 ) 


Based on the theory of [32], the probability that a molecule is released in time slot and 
absorbed in time slot is. 


A3-i+l)T 

p'{i,j)= g{r,t)dt 

J U-i)T 

= G(r, {j -i + l)r) - G(r, (j - i)t) 
Considering the ligand-based model, (12) can be revised as. 


p{iG) 


aQp'{iJ) 

b 


( 12 ) 


(13) 


4.2. Molecular classification 

With the probability given above, we could calculate the molecule quantity that are absorbed by the 
receiver in various slots. According to this character, the molecules absorbed during one slot could 
be divided into three portions, i.e., signal molecules, inter-symbol interference (ISl) molecules and 
noise molecules. Similar with traditional communication, we establish a relation formula for the 
three portions. 


P{14rx{^r) — — P{^sig(jtt: Clj.'j 4- ^ ^ 'Ur) T l^noiseijUr) — (14) 

In (14), the three portions of the right side respectively mean the expectative signal molecules, 
ISl molecules and noise molecules that are received during nf slot. We assume the corresponding 
quantities are respectively fcg, fci and ^ 2 , s.t., ^ 11=0 

4.2.1. Signal molecules The first term of (14) Nsig{nt,nr) indicates the signal molecules in the 
nl^ slot, which are transmitted in nf slot. The relation between rit and rir mainly depends 
on channel condition, transmission distance and slot length. While every solo molecule in the 
environment performs a random movement, the macroscopical distribution could be calculated, 
which is expressed by molecular absorbing probability (13). The quantity of the signal molecules 
follows the binomial distribution [10], 


piNsigint,nr) = ko) = ^ p(n4, (1 - p(n4, (15) 

Notes that (15) is the signal molecules of bit ”1”. If bit is ”0”, (15) becomes 0. 
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4.2.2. ISI molecules The inter-symbol interference (ISI) exists in the diffusive model also because 
of the randomness of molecular movement. As described before, each receiving slot maps a 
transmitting slot. We divide ISl into two themes, i.e., the interference from other bits and the 
interference to other bits. For a specific bit, 

• The first theme indicates that the molecules of the non-transmitting slots are absorbed by 
receiver in receiving slot, they are considered as the ISl from other bits. 

• The second theme indicates that the molecules of the transmitting slot are absorbed by receiver 
in non-receiving slots, they are considered as the ISl to other bits. 

The first theme increases the error rate of bit ”0” demodulation. Because some molecules of 
the other bits are added in the demodulation of bit ”0”, increase the probability of exceeding the 
threshold to decode it to ”1”. In a similar way, we could get that the second theme increases the 
error rate of bit ”1” demodulation. 

An example could be listed in Figure 2, the dash line from nt(2) to nr(l) is the first theme, 
indicates the ISl from second bit to first bit. The dash line from nt(l) to nr{2) is the second theme, 
indicates the ISl from first bit to second bit. On condition that the length of slots r is long, ISl of the 
adjacent bits is much heavier than that of the nonadjacent bits. So we only consider the ISl of the 
adjacent bits, which is presented by expanding the second item of (14), 


p{N{nt - 1, Tlr) + N{nt + 1, Ur) = fci) 

^ (16) 
= 2 _^P{N{nt - l,nr) = h)p{N{nt + l,nr} = h - h) 

h=0 

Flere, those of [nt — 1)*^ and (n* -I-1)*^ are regarded as the interference molecules as a view of 
time slot. To further expanded, the expression is. 


p{N{nt - 1, Ur) + N{nt + 1, Ur) = fci) 

h=o ^ ^ 

Note that ISl of Ijit -1-1)*^ bit exists only when nt <nr. The relation between nt and nr is 
determined by the transceiver, nt = Ur means that signal molecules which are emitted in one slot 
will be received in the same slot with a large probability, requiring that the distance between the 
transceiver is short or the slot length is large. 

4.2.3. Noise molecules The third item of (14) is the channel noise of the system. It is mainly caused 
by other nano-devices if nano-devices communication exist in the environment [28], apart from the 
possible molecular reaction and background molecules. When receiving and counting undesired 
molecules in demodulation, those molecules could be regarded as noise. In our literature, we cite 
the counting noise at the receiver proposed in [26]. The counting noise is accumulatively relevant 
to the quantity of nano-devices in the environment. More nano-devices will lead the heavier noise. 
On the contrary, if only one transceiver pair exists, the channel condition is much better. The direct 
impact of noise to the receiver is altering the quantity of the receiving molecules, i.e., the Nrx{i) in 
(7). Positive noise means receiving redundant molecules from other nano-devices. Negative noise 
means some signal molecules are received by other nano-devices. We assume that the noise is 
Additive White Gaussian Noise (AWGN), which is similar with [28] and expressed as, 

Nnoise{nr) ~ Normal{0, a^) (18) 


p{nt - 1, nr)'^{l - p{nt - 1, 


(17) 
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In (18), cr^ is the variance and indicate the noise power. It is positively related to the quantity of 
nano-devices and some other factors. Taking the expression of Normal distribution into account and 
the probability of third item in (14) is, 


noiseij^r') — ^ 2 ) 


1 

,- e 


(19) 


5. ATV ALGORITHM 


In this section, we first prove the feasibility of ATV algorithm in theory analysis, i.e., why we design 
ATV, based on the classified model of section 4. Then, we design the detail ATV algorithm based 
on the theory analysis, describe how it works. After that, we calculate the complexity of ATV, and 
discuss its performance on various channel conditions. 


5.1. Theory analysis 

The goal of the theory analysis is to prove the feasibility of ATV design. Theorem 1 below gives the 
theoretical optimal receiver threshold in the demodulation of diffusion-based MC. Considering 
the impact of the random noise, the Nip’* changes with the channel condition, we denote its mean 
as 

Proposition 1 

Let E{Nl^} and E{N^^} be the mean value of received molecules for bit ”1” and ”0”, then 

^ E{KJ+E{N°J 

Proof 

From (14) we could see that received molecules include signal, ISI and noise branches in the 
demodulation. Focusing on the noise term, Nnoise is a random Gaussian variable with zero-mean 
and variance of a^. Hence, for bit ”1”, we express the probability of absorbing Nrx molecules as, 

p{Nrx\l) = ^ e (20) 

sJ2'eu 

In a similar way, for bit ”0”, we have, 

p{Nrx\ 0 ) = r— e (21) 

s/2'eu 

Np’* should be the one that minimize BER performance in the demodulation, expressed by. 


NP* = argmin pe (22) 

N'j' £ [0, Af] 

So Nt = Np* requires that = 0, substituting (8) to (10), we have, 

aA'^Q 

pm = l)p{Np*\l) - pm = 0)p{Np*\0) = 0 (23) 

Substituting (20) and (21), we calculate Np* as. 


]Sl0 jyo jyl 

j^opt _ isi ' sig ' isi ' sig 

^ ~ 2 


_ p(6(i) = 0) 

NL + - Ng, - ''*pm = 1 ) 


-In 


(24) 


In typical case, the probabilities of emitting ”1” and ”0” are equal. And the mean value of Nnoise 
is 0, so we get its mean form as. 


^ E{NgJ + E{Ngp + E{NP} + E{Ngp EjN^} + L;{iV°J 
'^2 2 


(25) 

□ 
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Note that and are both variables that are determined by multiple factors, such as the 
channel noise, medium type, slot length, etc. So the is also a variable, that’s why we calculate 
its mean. Slot length is an important factor in the system, which is related with the ISI and delay. In 
this paper, considering its impact on ISI, we list the proposition 2 as follows. 

Proposition 2 

2 

If rit = n^, r > gjj is a needed not sufficient condition of N^ig > Nisi. 

Proof 

According to the definition of signal and ISI molecules in section 4, they both follow the binomial 
distribution. When rit = Ur, only the ISI from previous rit slots to later slots exists. So we get 
the mean value of Nisi as M{G{2t) — G(t)}, based on the character of the binomial distribution. 
Similarly the mean value of Nsig is MG{t). So, we have. 


Ns^g > N,s^ ^ 2G(t) > G(2r) (26) 

G(t) and G(2 t) could be got for specific r because it’s the error function or Gauss error function, 

2 

listed in (11). It’s easy to calculate that to satisfy (26), r > -^ should be required. □ 

We choose ^ applied in theorem 2 because it’s the time to peak of g{x,t), that is important 
in demodulation and slot length designing. In theorem 3, we will prove the feasibility of ATV 
algorithm, with the acid of proposition 1 and 2. 

Proposition 3 

NfF^ < ^, on condition that nt = n^. 

Proof 

Based on proposition 1, we get the expression of 7Vg^‘ as, 

_ ^{^lig + ^isi + ^sig + + ‘^^noise} 

Note that E{Nnoise} = 0, E{N^^g} = 0, calculating the distribution of the remaining three items 
based on the binomial distribution, we have. 


^ M{2G(2r) - G(r)} 


(28) 


2 

According to proposition 2, to make sure Nsig > Nisi, r > ^ is required. Under this condition , 
it’s easy to know that N^^* of (28) is smaller than 


□ 


5.2. A7Y algorithm design 

At the receiver, the threshold Nt is fixed generally for the demodulation of bit ”1” or ”0”. As section 
3 said, for typical demodulation, the molecules of one bit is received, if the quantity of receiving 
molecules N^x is greater than Nt, the estimation is ”1”. Conversely, the estimation is ”0”. 

A serious problem in MC is the communication reliability. Compared with traditional 
communication, biological channel condition of MC is more probabilistic, that may result in the high 
bit error rate (BER). The optimal threshold that enables the BER lowest, is different under 
the various biological channel conditions. Moreover, N!}?* is also changing with the various bits 
because of the noisy randomness character. The receiver doesn’t know how to set the best threshold 
at the beginning of the communication without knowledge of channel condition. So the half of 
the molecule emitted per bit M /2 is set as the threshold originally based on the typical experience, 
while it is not optimal that is proved in proposition 3. That’s the motivation to design ATV algorithm, 
which aims to improve BER performance by adjusting the threshold. Adopting ATV, the receiver 
threshold changes naturally based on knowledge of received bits, which are related with channel 
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Algorithm 1 Adaptive-threshold Variation 

1: Set iVT(l) = M/2, m = 0, no = 0, = 0, = 0 

2: for i=l,i^ time slot amount,i-i-i- do 

3: if Nrxi}) ^ Nrii) then 

4: receive a ”1” bit 

5: + Nrx{i), ni-l-t- 

6 : else 

7: receive a ”0” bit 

8 : + Nrx{i), no++ 

9: end if ^ 

10: Calculate A{i) = Nrii) — 

11: Calculate — Nxii) 

12: if A(i) — B(i) >/r then 

13: Nrii -b 1) = Noil) - 1 

14: else if A{i) — B{i) < —p then 

15: -/Vt(i- b 1) =-/Vo(*) + 1 

16: else 

17: iVT(i-b 1) =-VQ(i) 

18: end if 

19: end for 


condition. It is similar with the learning process, which makes receiver fit the realistic MC channel 
environment to demodulate properly. 

The ATV is shown in the table Algorithm 1. The threshold of slot is denoted by A^t(*)- It 
varies dynamically, according to the distance between the previous mean received molecules and 
the previous thresholds. First of all, we set the initial threshold A^t( 1) as M/2, and two counters are 
set up, respectively counting the quantity of bit 1 and 0. Under the equivalent probability of emitting 
1 and 0, the probability of receiving 0 and 1 is not equal, which will be shown in Figure 3. And the 
fixed threshold may result in the increase of the error probability, as the time lapses. Inspired by 
proposition 1, the optimal threshold should be the midpoint of the mean quantity of bit ”1” and ”0”. 
But they change with time, so we calculate current the distance between the current threshold and 
the two mean values of bit ’T””0”, expressed as A(j) and Bii). If the distance difference of A(i) 
and B{i) exceeds the tolerant interval p, we treat that the balance of bit 1 an 0 is broken and the 
threshold will accordingly be regulated. 

5.3. ATV algorithm discussion 

Imaged that nano-devices range from nanometer to micrometer scale. Restrained by the tiny size 
and low power, not many resources are available to implement the sensing, communication, acting 
functions. So to fit the character, any algorithms are required to be simple, efficient. ATV algorithm 
contains three major steps. The first step is to set the counter of bit ”1” and ”0”, and memory the 
molecule quantity per time slot, which ranges from line 1 to line 9. The time complexity of the first 
step is n -b 1 + nno + nni, which equals 2n^ -b n -b 1. The second step is to calculate the numerical 
distance between current threshold and average molecules received for bit ’T””0”, that ranges from 
line 10 to 11. The time complexity of the second step is 2n^. The third step is the threshold variation 
process, ranging from line 12 to 19, of which the time complexity is o?. Adding the three parts 
together, the entire complexity is bin? -b n -b 1. So We have the time complexity of ATV as. 

Tin) = Oin^) (29) 

In ATV algorithm, only one ’’for loop” is used. Compared with some other existing methods 
which improve reliability in molecular communication, the best advantage of ATV is its 
simplification, low time complexity, that make it more likely to be utilized in nanoscale devices. 
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Table 1. SIMULATION PARAMETERS 


Diffusive Coefficient 

D 

10,1000 

j s 

Transceiver distance 

r 

1-20 

p 

Time slot length 

r 

1-10 

s 

Binding rate of the receptors 

a 

0.1 

1 

Releasing rate of the receptors 

b 

0.08 

1 

Concentration of the receptors 

Q 

1 

umol/I 

Power of AWCN 

p . 
noise 

1-20 

pW 

Receiver threshold 

Nx 

50-450 

/ 

Molecules for bit ”1” 

M 

500 

/ 

Tolerant interval 

M 

30,60 

/ 


We can infer that the varying threshold is convergent with a long enough time. The reason is that the 
signal and ISI components of N^x approach their fixed means when calculating and iV°^. The 
random impact of the noise term to A(i) and B{i) becomes smaller and smaller with the increase of 
the Nxx- So A(i) and B{i) have a limit, that makes the variation of the threshold convergent. 

ATV algorithm aims to keep the balance of demodulation, which are able to reduce the error 
probability of demodulate bit ”1” or ”0”. Channel noise is an important factor to rise the BER. As 
we describe before, the noise is assumed as the AWGN for demodulation. Serious noisy channel 
here means a noise with a larger variance, which may be caused by crowded communication nano¬ 
devices. Mapping into the mathematical models of (14), the receiver is possible to absorb or lose 
more molecules. Hence the value of A{i) — B{i) m ATV has a larger variable range, which makes 
Nxii) easier to change. Stronger noise also means a lager random factor, that lead the delay of 
threshold converge. Differently, in a good channel condition, there only exists slight noise. If the 
slot length is chosen properly, compared with signal, ISI is relatively weaker in terms of power. The 
threshold variation will not be fiercely and is easier to converge. 


6. PERFORMANCE EVOLUTION 

In this section, we simulate the reliability of free diffusion-based channel detailed in previous 
Sections. Considering the ISI and channel noise, BER of the channel is presented with some 
important parameters. The performance of the ATV algorithm is also verified. Note that we only 
consider the case that n* = Ur. Let 7 e be the signal to interference ratio (SINR), that is calculated 
based on the molecule ratio of the three branches in Section 4, expressed as. 




7e = 


signal 


^El=iW{nt,nr)\‘^ 


^ El=i + 1) rirW + 


(30) 


Where Psignai, Pisi and Pnoise respectively means the power of received signal, ISI signal and 
channel noise. They are established like this because molecule quantity is a discrete value. [.] is 
the floor operator. (15) reflects the ratio of each part of molecules directly. When receiver counts 
molecules in every time slot, ISI and channel noise molecules result in the quantity inaccuracy 
and promote error bit determine. And they could be considered as the addictive signals comparing 
with receiving signals according to the character of concentration modulation and quantity threshold 
decision, n is the total quantity of the time slot, and Pnoise is its variance since it follows a zero-mean 
Normal distribution. 
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{a} r = 4|im, D = 1 0 jim/s (b) r = 4)im, D = 1 000 jim/s (c) r = 20|im, D = 1 000 ixm/s 





Figure 3. BER of communication for different D, r and r under various 7 e 



Figure 4. BER and Error rate of bit r=4s, D=10/rm^/s, r=4/rm, 7e=10 


6.1. Simulation Design 

The major simulation parameters are listed in table 1. These parameters could be divided into two 
branches, i.e., diffusion-based channel parameters and signal (including interference and noise) 
parameters. In [31], the diffusion coefficient D is considered as the range of 10-1000 [irr? js. The 
channel conditions under different D are verified. In this paper, we study the diffusive channel based 
on some discrete D. We set the transceiver distance r as 1-20 gLm, similar with the references [29] 
and [31]. The value of a,b and Q indicate the parameter of the ligand-based model, that are derived 
from [21]. In our paper, molecules quantity per bit M is set as 500. The setup of /r mainly refers to 
the result of Figure 5 and 6. We set it as 30 and 60, when the time slot length is 2s. 


6.2. Numerical analysis 

Figure 3 illustrates the relationship between BER pe and channel diffusion D, transmission distance 
r, as well as time slot length r. The varies through altering the noise power. We can see 
that Pe decrease with 7 e increase. There is little difference for pe of different r when 7 e is low. 
Comparing (a) and (b), we could get the conclusion that the BER performance does not have a 
significant difference for different D. While comparing (b) and (c), larger the distance r is, worse 
the BER performance is. This attributes to two factors. Eirstly, increase of r directly results in 
decline of signal molecule arriving probability, which refers to (3). Correspondingly, the impact to 
ISl molecule is no far more significant than signal molecules. Moreover, for a longer r, the signal 
molecules decrease obviously, so channel noise will impact more when receiver counts molecule 
per time slot. 
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Figure 5. BER for various and r, D=\Q^im?/s, r=4fj,m, 7e=10 



Figure 6. BER for various Nt and r, /s, r=4s, 7e=10 


Based on Figure 4, we can see that and p\ decrease and increase respectively with the increase 
of Nt- The crossing point of the two curves decides the optimal value of Nq. It is smaller than the 
midpoint, i.e., according to the simulation. If current iVo is much greater than the optimal 

value, Pg falls lightly while p\ rises heavily, as a result the BER becomes large. We can get the 
similar conclusion if Nt is much less than the optimal value. 

In Figure 5, the BER pe is shown with various Nt under different time slot r. pe first decreases 
and then increases, with enlargement of the threshold. For different T, the optimal threshold N!^^ 
is different. The molecule quantity per bit M is set as 500, while N^^ is lower than half of the 
molecule amount per bit, i.e. 250, owing to the specificity of molecular OOK modulation. It also 
presents that N^^ trends to get a greater value with the increase of r. Correspondingly, Figure 
6 reveals the relationship between BER and threshold under different r. The variation trend for 
threshold is the same with Figure 5. However, Pf. is lower for lager r, and N!p‘* is also farther from 
the middle point. Hence, we can conclude that for various r and r, the threshold for minimum pe is 
various. 

In Figure 7, the BER performance of using or not algorithm 1 under various r is demonstrated. 
Although when 7 e is extremely low, the performance improvement is not obvious, even worse. 
However, when 7 e is higher, use of ATV algorithm could lead a distinct difference. The receiver 
could not obtain the right amount of the signal molecules and the proper threshold will be deviated 
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Figure 7. Comparison of using or not using ATV under various 7 e, D=lOfj,m/s, r=Sfj,m 


(a) (i=30 




Figure 8. Receiver threshold varied with time, r=2s, D=lOfim^/s, 7 e = 10 


to the original value. The ATV algorithm makes the actual threshold closer to theoretic optimal 
threshold, so the pe decrease. In addition, for a shorter r, adopting ATV could get a more significant 
BER performance, gaining smaller channel delay. 

Figure 8 and 9 show the variation of receiver threshold if the ATV algorithm is adopted. In those 
figures, we can see that the threshold decreases sharply firstly. There exists much sawteeth due to 
the randomness of AWGN and molecular Brownian motion. Comparing figure 6 (a) and (b), the first 
curve shakes a little heavier, just because that the tolerant interval p is smaller, leading threshold 
changing easier. For Figure 8 (a) and Figure 9 (c), only the condition of 7 e is different. We could 
see that, lower 7 e, i.e., worse communication environment, causes much more shake of the varied 
threshold. So the pe will also increase accordingly. We could conclude that in both good and a little 
worse communication condition, the threshold of ATV is more closer to best theoretic threshold. 
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(c) H = 30 




Figure 9. Receiver threshold varied with time, r=2s, D=l0fim^/s, 7 e = 5 


7. CONCLUSIONS 

In this paper, we investigate the reliability of the diffusion-based molecular channel. We first 
describe the basic system model. A classified model is established to divide molecules into signal, 
ISI and noise parts in demodulation. We expand the three parts respectively and connect them using 
a traditional-like formula. Based on the classified model, we do a theoretical analysis, which give the 
basis for the later ATV algorithm design. The ATV algorithm at the receiver is designed to improve 
the BER performance. We discuss its complexity and performance on various noisy channels. Then, 
we verify our proposal through simulation. The results show that the algorithm decreases BER 
obviously as long as the channel condition is not too bad. The threshold of the receiver is more 
closer to the best theoretic threshold, so the algorithm works. This work contributes to improve the 
communication reliability for diffusion-based molecular communication. 
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